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A  framework  for  evaluating  the  epitaxy  of  ordered  organic  overlayers  of  generic  symmetry  on  ordered  substrates  is 
described  that  combines  a  computationally  efficient  method  for  explicit  determination  of  the  type  of  epitaxy  (i.e., 
commensurism,  coincidence,  or  incommensurism)  and  overlayer  azimuthal  orientation  with  an  analysis  of  the  elastic 
properties  of  the  overlayer  and  the  overlayer-substrate  interface.  The  azimuthal  orientations  predicted  by  this  function 
agree  with  values  predicted  by  semi-empirical  potential  energy  calculations  and  observed  experimentally  for 
previously  reported  organic  overlayers  which  are  demonstrated  here  to  be  coincident,  including  electrochemically 
grown  overlayers  of  molecular  conductors.  Calculations  based  on  this  analytical  approach  are  much  less 
computationally  intensive  than  potential  energy  calculations  as  the  number  of  computational  operations  is 
independent  of  the  overlayer  size  chosen  for  analysis.  This  enables  analyses  to  be  performed  for  the  large  overlayer 
basis  sets  common  for  naolecular  overlayers.  Furthermore,  this  facilitates  the  analysis  of  coincident  overlayers,  for 
which  the  overlayer  size  needs  to  be  large  enough  to  establish  the  phasing  relationship  between  the  substrate  and 
large  non-primitie  overlays  supercell  so  fi\at  the  global  minimum  with  respect  to  azimuthal  angle  can  be  determined. 
The  computational  efficiency  of  this  method  also  enables  convenient  examination  of  numerous  possible  reconstructed 
overlayer  configurations  in  which  the  lattice  parameters  are  bracketed  around  those  of  the  native  overlayer,  thereby 
allowing  examination  of  possible  epitaxy-driven  overlayer  reconstructions.  When  combined  with  calculated  intralayer 
and  overlayer-substrate  elastic  constants  this  method  provides  a  strategy  for  the  design  of  heteroepitaxial  molecular 
films. 
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version  2.0:  May  3,  1995  between  the  overlayer  and  substrate.*  Controlling  these  factors  is  imperative  as  reconstruction  of 

tCiuient  address:  Depjuunent  of  Chemical  Engineering.  University  of  Virginia.  Charlottesville. 

VA  22902.  the  primary  overlayer  from  its  native  form  or  stress-induced  defects  can  affect  the  quality  of 

PACS  61.16.Ch;  62.20.Dc;  68.55.-a;  81.15.Lm  multilayer  Alms  and  bulk  crystals  grown  from  the  primary  overlayer. 


While  there  have  been  significant  advances  in  the  understanding  of  the  physical  and  fornt  The  ntetiiod  also  enables  convenient  searching  for  epitaxial  relationships  between  a  rigid 

electronic  properties  of  molecular  films,  lliere  exists  a  need  for  paradigms  enabling  a  priori  design  substrate  and  many  possible  reconstructed  forms  of  an  ovcrlayer,  aiid  suggests  liiat  a  priori  design 

of  heteroepitaxial  molecular  overlayers.  We  have  been  en^loying  a  strategy  for  the  design  of  of  overlayer-substrate  systems  is  feasible, 

molecular  films  in  which  the  native  structure  of  a  molecular  overlayer,  considered  to  be  at  or  near 

Its  minimum  energy  configuration,  is  surmised  from  iu  structure  in  bulk  crystals,  which  generally  METHODS  AND  MATERIALS 
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coincident  ovetUyers  for  which  reeUtry  with  substrates  can  only  be  established  by  large  non-  containing  the  respecuve  eleettolytes  under  condiuons  similar  lo  those  previously  repotted  for  the 

ptimiuvesupetcells.  When  combined  with  calculations  of  overlayer  and  overlayer-substrateelasUc  electrosynthesis  of  bulk  ctystab  S*-'*-**  Au(lll)  subslrates  for  deposiUon  of  the  (TTF)(TCNQ) 

constanu  for  native  ovetlayets  whose  structures  are  surmised  from  either  crystal  stiuctutes.  ovetlayen  were  prepared  by  melting  high  purity  gold  wire  (99  999%)  in  a  oxygendiydtogen  flame 

calculations  or  experimental  data,  this  approach  enables  rapid  analysis  of  overlayer-substrate  to  expose  (1 11)  oriented  facets, 

systems  and  a  qualititative  assessment  of  the  tendency  for  overlayer  reconstruction  from  its  native 


Scanning  tunneling  (STM)  and  atomic  force  microscope  (AFM)  experiments  were  personal  computer  using  a  program  written  in  our  laboratory  (EpiCalc®),  which  runs  in  U»e 

performed  with  a  Nanoscope  III  Multimode  scanning  probe  microscope  (Digital  Iristruments,  Santa  Windows^^  v.  3.1  environment.  EpiCalc®  is  available  on  the  World  Wide  Web  at 

Barbara,  CA).  STM  tips  consisted  of  mechanically  cut  Pl/lr  wires  and  AFM  probes  (Nanoprobe,  hltp://www.cems.unm.edu/research/ward. 

Park  Scientific,  Palo  Alto.  CA)  consisted  of  triangular  silicon  niuide  cantilevers  (force  constant  «• 
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The  primary  value  of  the  iraiisformalion  matrix  is  tlial  it  describes  the  interface  in 
The  Epitaxial  Interface.  Epitaxy  generally  is  used  to  describe  lattice  registry,  or  equivalently  convenient  terms.  The  determinant  of  C  is  equivalent  to  Uie  ratio  of  overlayer  to  substrate  unit  cell 

me  degree  of  "phase  matching",  between  two  opposing  lattice  planes  (aimough  it  was  originally  areas  and  its  value  defines  whemer  the  system  is  commensurate,  coincident,  or  incommensurate, 

used  to  describe  the  growth  of  material  B  on  subsUate  A,  with  B  adopting  me  structure  of  A).  In  However,  mis  determinant  does  not  provide  the  energetics  of  die  interface.  Because  the 
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commensurate  wim  the  substrate.  This  is  equivalent  to  statmg  that  a  comcident  ovcrtayer  nas  some 
overlayer  positions  (most  commonly  chosen  to  be  me  vertices  of  me  supercell)  which  are 
conanensurate  wim  substrate  sites,  while  omer  molecules  contained  within  me  supercell  are  locally 
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one-diniensional  lattice  can  be  described  by  eq.  (9),  where  a  represents  the  lattice  periodicity. 


llie  potential  WNo  can  be  calculated  for  orientations  defined  by  the  component  of  maU^ix  C  indicates  that  the  dependence  of  Vt  on  6  is  exactly  reproduced  with  VA^o  for  this  overlayer> 

for  an  overlayer  with  M  x  N  unit  cells.  This  potential  assumes  discrete  values  which  depend  upon  substrate  interface.  Actual  molecular  overlayers  are  likely  to  exhibit  differences  in  Vt  and  V/Vq  as 

the  type  of  epitaxy,  with  commensurism.  coincidence,  and  incommensurism  resulting  in  values  of  VA/  o  does  not  account  for  local  molecule>substrate  interactions.  However,  as  the  over  layer  size  is 
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The  total  potential  Vt,  two-dimensional  misfit  D(d)i  VA^o  for  a  trivial  5x5  overlayer  calculated  from  V/Vq  using  the  known  lattice  parameters  ffable  1).  In  nearly  every  example  in 

and  substrate,  in  which  aj  =  fl2  =  bj  =  62  and  a  =  p  =  90®,  each  exhibit  minima  at  orientations  of  6  Table  1.  the  orienUlion  of  the  overlayer  prediaed  by  V/Vo  corresponds  to  that  determined 

=:  0®  and  0  =  90®  (Figure  3)  conesponding  to  a  commensurate  condition.  The  shape  of  D(®)  experimentally.  Particularly  interesting  are  the  six  polymorphs  of  perylene-tetracarboxylic 


snail.  The  emergence  of  a  clear  global  minimum  associated  witli  coincident  overlayers  for  large  interactions  will  be  maximized  by  formation  of  a  conmiensuraie  overlayer,  or  to  a  lesser  extent  by 

overlayer  sizes  denansuates  that  multiple  minima  or  shallow  potential  hmctions  observed  for  small  coincidence,  as  Uie  phasing  of  the  potential  functions  is  optimized  under  these  conditions, 

overlayer  sizes  does  not  necessarily  signify  incommensurisia  The  concept  of  •  quasiepitaxy,*’  a  Intralayer  forces  favor  retention  of  the  native  form  of  Uie  overlayer  wlule  overlayer-substrate  forces 

condition  surmised  from  apparent  lack  of  comii^nsurisin,  was  advanced  recently  to  explain  the  at  the  interface  promote  reconstruction  of  the  overlayer  m  order  to  attain  commensurism  or 


greater  dispersive  interaction  per  molecule.  Consequently,  the  interfacial  elastics  constants  are 
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n¥)l“'  (at  e  =  0®)  and  the  variation  in  thia  energy  is  0.8  kcal  mol*',  both  substantially  higher  than  with  such  an  overlayer  is  not  possible,  its  large  intralayer  elastic  constant  will  uihibit  the 

that  calculated  for  the  (OOi)  p-(ET)2l3  overlayct.  This  is  due  to  sttonger  interfacial  interactions  reconstruction  required  to  achieve  coincidence  or  commensurUnt  Under  this  condiuon  the  system 

associated  with  the  Pe  molecular  planes  lying  parallel  to  the  HOPG  surface,  which  results  in  a  will  likely  favor  the  growth  of  epitaxial  overlayers  whose  structures  resemble  alternative  bulk 


crystal  planes  with  smaller  intralayer  interactions  and  elastic  constants.  This  reduces  tlie  energetic  overlayers  serve  as  nuclei  for  multilayers  or  bulk  crystals,  tiiis  can  serve  as  a  convenient  approach 

penalty  of  overlayer  reconstruction  necessary  to  achieve  epitaxy  by  either  coincidence  or  to  selecting  substrates  for  selective  crystallization  of  tlnn  film  structures  or  crystals  that  lend  to 

commensurism.  However,  this  can  lead  to  growth  orientations  in  which  the  suongest  exhibit  polymorphism.  Wlien  combined  with  analysis  of  the  elastic  constanU  for  a  given 

intermolecular  interactions  are  perpendicular  to  the  overiayer,  as  is  the  case  for  (Pe)2C104llH0PG,  overlay er-substrale  configuration  this  provides  a  convenient  approach  for  the  a  prior  design  of 

(Oil)  O-fEDihliHOPG,  and  (TTF)(TCNQ)liAu(Ul).  Under  this  condition,  the  elastic  constants  heteroepitaxial  ntolecular  films  and  substrates  for  the  growtli  of  multilayer  films  and  bulk  oystab. 
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used  to  determine  the  likelihood  of  epUaxy  of  different  polymoriduc  forms  of  an  overiayer  if  these  ^  (a)  G.E.  Collins,  K.W.  Nebesny,  C.  England,  L.  Chau,  P.A.  Lee,  B.  Parkinson.  Q. 

polymorph  structures  are  known,  providing  son»  predictability  of  the  overlayer  structure.  If  these  Fernando,  and  N.R.  Armstrong,  J.  Vac,  Sci.  Technol.  A  10,  2902  (1992).  (b)  C. 


Ludwig,  B.  Gompf,  W.  Flatz.  J.  Peterson.  W.  Eisenmenger.  M.  Mobus.  U.  20  a. I.  Kitaigorodsky.  Molecular  Crystals  and  Molecules,  (Academic  Press.  New  York, 

Ziinmermann,  and  N.  Karl,  Z.  Phys.  B  86,  397  (1992),  (c)  C.  Ludwig,  B.  Gompf,  J.  1973). 

Peterson,  R.  Strohmaier,  and  W.  EisenirKngcr.  Z.  Phys.  B  93.  365  (1993).  (d)  N.R.  21  MM2/MMP2:  (a)  N.L.  Allinger,  J.  Am.  Chem.  Soc.  99.  8127  (1977).  (b)  J.T.  Sprague, 
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Figure  Captions  (c)  (Pe)2CI04llH0PG,  (d)  (TrF)(TCNQ)llAu(l  11).  (e)  P-(ET)2i3  (Type  1)11  HOPG. 

Figure  1.  Schematic  representalion  of  a  generic  2x2  molecular  overlayer  on  a  rigid  substrate.  p-(ET)2l3  (Type  11)11  HOPG.  The  calculations  were  perfornted  using  20  x  20 

The  substrate  and  overlayer  lattices  are  defined  by  two-dimensional  cells  with  lattice  overlayers. 

constants  aj,  a2,  a  and  62,  p,  respectively.  The  angle  8  represents  the  angle  Figure  6.  Schematic  representations  of  overlayer  orientations  on  substrates  for  selected  systems 
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1  atoms  per  unit  cell.  was  fixed  at  the  mimimum  established  in  (b).  Tlie  local  minimum  at  6  ~  19^  reflects 

Figure  5.  Calculated  dependence  of  V/Vq  on  azimuthal  angle  0  for  the  selected  overlayer-  ^  minimum  observed  for  large  overlayer  sizes  (see  Figure  4),  whereas  the  deeper 

substrate  combinations  depicted  in  Hgure  6.  (a)  PrCDA)ICu(lOO).  (b)  CuPcliMoS2*  minimum  evident  at  60^  is  a  false  mininuim  resulting  from  the  small  overlayer  size 


used  in  the  calculation.  (e.Q  Dependence  of  V,  Ox.  Oy.  Cxx  and  Cyy  for  for  an  isolated 
the  (001)  p-(E02l3  layer.  The  minima  of  the  potential  wells  are  a  =  6.67  A  and  a  - 
9.09  A,  essentially  identical  to  the  conesponding  values  of  a  =  6.61  A  and  6  =  9.1  A 
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